Efficiency scaling of non-coherent upconversion 
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A very promising approach to obtain efficient upconversion of light is the use of triplet-triplet annihilation 
of excitations in molecular systems. In real materials, besides upconversion, many other physical processes 
take place - fluorescence, non-radiative decay, annihilation, diffusion - and compete with upconversion. The 
main objective of this work is to design a proof of principle model that can be used to shed light on the 
relevance of the interaction between the different physical processes that take part in these kinds of systems. 
Ultimately, we want to establish general principles that may guide experimentalists toward the design of 
materials with maximum efficiency. Here wc show, in a ID model system, that even in the presence of these 
processes upconversion can be optimized by varying the ratio between the two molecular species present in 
this kind of materials. We derive scaling laws for this ratio and for the maximum efficiency of upconversion, 
as a function of the diffusion rate J, as well as of the creation and of the decay rate of the excitations. 



I. INTRODUCTION 

To obtain highly efficient solar cells is one of the most 
important technological problems to be solved in the fu- 
ture. In the last years, researchers have developed new 
strategies to achieve this task: solar cells from novel 
materials, multi-junction solar cells, but also, new and 
cheaper fabrication techniques, etcPES 

Recently, to overcome the efficiency limit due to spec- 
tral mismatch^, the idea of upconverting low energy pho- 
tons has been revivecP^. The idea is to absorb those 
photons that can not be absorbed by the solar cell, and 
send them back to the cell with an appropriate wave- 
length, ie. a wave length that can efficiently be absorbed 
by the cell. One example of a well studied low power 
upconversion mechanism is noncoherent upconversion by 
triplet-triplet-annihilation (TTA^SH^H. The migration of 
excitation energy through assemblies or aggregates of or- 
ganic light absorbing molecules is used for capturing the 
energy of sunlight, similar to organic solar cells and nat- 
ural light harvesting systems 2 ^", and is then transferred 
to a molecule where the upconversion occurs. The high 
absorptivity of these organic molecules, ie. in compari- 
son to nanoparticles^, is beneficial for the upconversion 
process. Other techniques, such as high harmonics gener- 
ation or excited state absorption, require high light inten- 
sitie s and usually exhibit only low upconversion quantum 
yielcpEIl. 

Successful experimental implementations of non- 
coherent upconversion exist for various materials^" 2 ^, in 
solutions and in thin films. To the best of our knowledge, 
the highest upconversion yield achieved so far is 19.3%, 
reported by Ji et alP^*. However, only few studies investi- 
gating concentration and power dependence exist! 13 * 18 * 26 *, 



and little is known about the quantitative effect of the 
various loss channels on the upconversion yield. In the 
present report we study a simple model system in order to 
address the question how the upconversion efficiency can 
be optimized by varying the ratio between sensitizers to 
acceptors - the two molecular species that are present in 
materials used for upconversion by TTA (see below) . We 
carried out model calculations to investigate the effect of 
sensitizer to emitter ratio, triplet lifetime and triplet cre- 
ation rate on the upconversion efficiency. We find that 
for the optimal choice of sensitizer-to-emitter ratio the 
efficiency can reach almost unity, neglecting spin statisti- 
cal effects. This leads us to suggest that optimal systems 
for non-coherent upconversion can be very efficient, when 
tailored to the light conditions of their environment. The 
results of our model give insight into dynamics of a non- 
coherent upconversion material and extend our intuitive 
understanding of the competition between different reac- 
tion pathways. 

The state of the art scheme f or up conversio n by TTA 



involves two types of molecule d 11 * 12 *, see also Figure 1 
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The sensitizer molecules harvest solar energy. Upon pho- 
ton absorption a sensitizer molecule is excited from the 
ground state So to its first excited singlet state, Si- This 
is followed by rapid intersystem crossing (ISC^ to a 
triplet state, T\. The T\-Sq transition is spin forbidden 
and therefore suppressed, but the transition can be medi- 
ated by spin-orbit interaction. Hence, these triplets have 
a relatively long lifetime, which enables a large triplet 
diffusion length. Eventually, the triplet is transferred to 
the second moiety, the emitter molecules (sometimes also 
referred to as acceptors). The emitter molecule has the 
property that the level spacing between S\ and ground 
state equals approximately twice the spacing of its triplet 
and ground state. Therefore, it is energetically possible 
to pool the energy of two triplets in the emitter Si state. 
Thus, if an emitter is in its Ti state and another triplet, 
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from a sensitizer or emitter, is close by, the emitter's Si 
state can be populated. If both participating molecules 
are of the same s pecies , this is referred to as homo, else 
as hetero process |22129| A strong Si-Sq coupling conse- 
quently leads to emission of an upconverted photon. In 
fact, fluoresce nce from the emitter Si state is experimen- 
tally observed^™. 

However, other processes compete with this upconver- 
sion mechanism: A triplet-triplet encounter can not lead 
to upconversion if both of the two molecules are sensitiz- 
ers, since the sensitizer Si state does not have the appro- 
priate energy. Instead, if both molecule s are sensitizers, 
the well known triplet-triplet-annihilatiorP2^J T+T — > T 
occurs, which in statistical physics literature is referred 
to as coagulation and has been studied extensively^^. 
Consequently, a triplet-triplet encounter can only lead to 
upconversion if at least one of the involved molecules is 
an emitter. If that is the case, and the spin states of the 
two triplets add to 7715 = 0, upconversion can take place: 
T + T -*hv. 

Transport of triplets in the material can occur either 
by diffusion of the molecule itself, in case of solution 
based materials^, or by triplet hopping between neigh- 
bouring molecules. In our model system, we consider a 
solid state material, where transport of triplets occurs 
exclusively by triplet hopping which can be described by 
Miller- Abrahams rate d 37 * 3 ^ , see section II below. In total, 



we thus need to consider five processes to capture the full 
kinetics: Triplet creation in sensitizer molecules, triplet 
hopping between neighbouring molecules, coagulation (if 
two triplets meet at a sensitizer molecule), upconversion 
(if two triplets meet at an emitter molecule), and triplet 
decay via phosphorescence. In principle, it is straightfor- 
ward to write down a master equation for these processes, 
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FIG. 1. Sketch of level structure for (a) sensitizer and (b) 
emitter molecules. Photonic transitions are indicated by solid 
lines. Dashed lines represent transitions induced by triplet 
energy transfer, (a) Upon photon absorption, the excitation 
energy fku is transferred from the singlet Si state to the triplet 
Ti state, via intersystem crossing (ISC). This process occurs 
with almost unity quantum yield in molecules used for non- 
coherent upconversion. (b) The emitter level structure allows 
for excited singlet state occupation by pooling the energy of 
two triplets. Thus, the emitter molecule can emit a photon 
of energy hu' > fio; via fluorescence. 



however, solving it is a daunting task. 



Similar systems were studied before, for example on- 
site annihilation and fu sion i n one-dimensional disor- 
dered media was studied 221 ^ 1 by mapping the master 
equation of the corresponding classical random walk onto 
a quantum Hamiltonian. Schiitz et alP^l found an approx- 
imation to the ensemble averaged time evolution of the 
many body random walk, by studying the many parti- 
cle random walk in its dual environment, where calcu- 
lations simplify for specific Hamiltonians. In principle, 
it is also possible to implement coagulation, particle de- 
cay and creation in the Hamiltonian. Unfortunately, the 
problem no longer simplifies in the dual environment if 
the Hamiltonian contains particle creation. 
Ben- Avraham et alP^ studied particle diffusion on a one- 
dimensional, translationally invariant chain including co- 
agulation and particle creation. Their work is based on 
the method of intervals or inter particle distribution func- 
tion. This function gives the probability that a subunit 
of the system contains no particle, as function of the 
size of the subsystem. The benefit of this framework 
is that it gives a good intuition of the system dynamics. 
The latter can be described by a differential equation 
in the inter particle distribution function, which can be 
solved using standard methods^. In infinite systems, 
this method gives an exact solution. However, applica- 
tion of this method for other boundary conditions is not 
trivial. In particular, it is not possible to use it for finite 
systems with absorbing boundaries and particle creation: 
The mean density vanishes on the absorbing boundary, 
but is finite inside the system, due to particle creation. 
Hence, the mean density clearly must be a function of the 
spatial coordinate, and, consequently, the system is not 
translationally invariant and the method can not be ap- 
plied. One important result - suggested already earlier^, 
but confirmed by Ben- Avraham et alP^ - is that coagu- 
lation in one-dimensional systems scales with the mean 
density cubed. We will use this fact below. 



Because of these difficulties, we will follow a more in- 
tuitive approach, looking for analytical solutions only in 
the regime in which the equations involved are linear. 
We will then deduce general scaling arguments, for the 
linear and nonlinear regime, from first principles. All our 
results will be underpinned with numerical simulations. 



The remaining four sections of this paper are struc- 
tured as follows: In the next section we introduce the 
upconversion model system. Then, we present numerical 
findings for the upconversion efficiency in this model. Af- 
terwards, we perform an analytical analysis for the mean 
field density and the upconversion efficiency. Finally, we 
present an intuitive argument that results in scaling laws 
for a critical system parameter, and for the upconversion 
efficiency. 
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II. THE MODEL 

To enhance energy transfer from sensitizer to emitter, 
one should choose an emitter with a triplet energy level 
below that of the sensitizer triplet state. In this case, 
energy transfer from emitter to sensitizer is unlikely to 
happen, since the corresponding Abrahams-Miller tran- 
sition rate is exponentially suppressed for AE ^> k B T, 
with temperature T and energ y difference AE > be- 
tween the respective level$*Z122l. This is the case for con- 
figurations where high efficiency was observed: In the ex- 
periment reported by Cheng et al™ rubrene was used as 
emitter, with a lowest trip let st ate energy of m 1.15eV, 
and a lifetime of sa lOO/iiPM The sensitizer PdPQ 4 
was excited with 670nm laser pulses, corresponding to 
fa 1.8 eV, ie. much higher than the energy of the T\ 
emitter state. Since AE = (1.8 - 1.15) eV = 0.65 eV > 
ksT ~ 0.03 eV for room temperature, it is reasonable 
to neglect emitter to sensitizer energy transfer in our 
model. On the other hand, we describe sensitizer to 
emitter transfer and transfer between two sensitizers by 
a constant rate J. This rate is determined by the inter- 
molecular distance and the electronic coupling between 
the respective states, which we assume to be identical 
for all neighbouring molecules. The lifetime of emitter 
triplets, « 100/xs, is, when compared with other relevant 
experimental values used by Cheng et alP^, the slowest 
timescale in our model. Hence, it is well justified to as- 
sume an infinite lifetime of emitter triplets in our model. 
As further simplification, the diffusion limited reactions 
upconversiorPS and coagulation are modelled as on-site 
interaction with infinite rate. Thus, two triplets occu- 
pying the same site immediately undergo upconversion 
(coagulation), if the site is an emitter (sensitizer), see 
below. Finally, the lifetime of the sensitizer triplet de- 
pends strongly on the material, values from 0.25ns up to 
lOO^s were reported^. 

As a consequence of the simplifications introduced 
above, excitations are trapped at emitters. Thus, in our 
one-dimensional model we can think of emitters as ab- 
sorbing barriers. This constraint simplifies the kinetics. 
As we will see in |section IV} it allows us to reduce the 
study of dynamics to sensitizers only, since each triplet 
which escapes to an emitter will eventually contribute 
to emission of an upconverted photon. In our model, we 
consider only hetero type upconversion, because this type 
yields higher efficiencies for otherwise identical parame- 
ters, and is therefore more promising. 

Our one-dimensional toy model is illustrated in the 
following sketch, where E and S represent emitter and 
sensitizer molecules, respectively. 
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The system of length L is characterized by the inter- 
molecular spacing Ax and the emitter spacing I Ax. 



Thus, the parameter I controls the sensitizer-to-emitter 
ratio. The system can be grouped into n blocks, each 
containing only a single emitter and I sensitizers. 

In an ensemble of emitter and sensitizer molecules, 
the energy of the Ti state will differ from molecule to 
molecule, due to dynamical disorder. This on-site energy 
disorder can be neglected, since for symmetric hopping 
rates it merely corresp onds t o a rescaling of the diffusion 
constant D = {Axf jinHEHl l n the time average, the 
on-site energy differences vanish and sensitizer hopping 
rates can be considered symmetric. Thus, we do not 
expect a considerable change in our results due to this 
approximation. However, this rescaling will decrease the 
diffusion constant and, thus, the range of validity of the 
low density approximation decreases as well. 

We perform an exact simulation of the upconversion 
system as depicted in ([I]), using the aforementioned ap- 
proximations and a kinetic Monte Carlo algorithm™^. 
|Figure 2| lists the processes included in the simulation: 
(a) triplet creation by absorption of a photon by a sen- 
sitizer molecule, with rate T c - note that an excited sen- 
sitizer does not absorb photons; (b) sensitizer triplet de- 
cay, with rate T^, e.g. via phosphorescence; (c) triplet 
hopping, with rate J between two neighbouring sensitiz- 
ers; (d) triplet energy transfer from sensitizer to emitter, 
modelled with the same rate J as triplet hopping among 
sensitizers; (e) triplet energy transfer, with rate J, im- 
mediately followed by coagulation/triplet-triplct annihi- 
lation; (f ) triplet energy transfer, with rate J, resulting in 
population of the emitter Si state, immediately followed 
by emission of an upconverted photon. 
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FIG. 2. Possible processes occuring in our upconversion 
model system. Black bars indicate the level structure, black 
dots represent occupied states, whereas arrows indicate tran- 
sitions. Sensitizers are represented by the two levels So and 
Ti, emitters by the three levels So, Ti and Si. The processes 
and respective rates are (a) triplet creation, with rate r c , (b) 
triplet decay, with rate Td, (c) triplet hopping, with rate J, 
(d) triplet energy transfer from sensitizer to emitter, with rate 
J, (e) triplet hopping with subsequent triplet-triplet annihi- 
lation or coagulation, with rate J, (f) triplet hopping with 
subsequent upconversion (not indicated), with rate J. 
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To measure the efficiency of the upconversion process, 
the internal quantum efficiency is a natural candidate: 



Xiqe 



of upconverted particles 
o. of created particles 



(2) 



However, this is a bad measure for the overall efficiency of 
a given configuration. It is only limited by the internal 
loss channels decay and coagulation. The efficiency of 
the absorption process is not taken into account. For 
example, consider a system which consists of only a single 
absorbing site, and of traps (emitters) otherwise. This 
system has Xiqe = 1, but will absorb only a tiny fraction 
of the available incident intensity. 

On the other hand, the external quantum efficiency, 
EQE, also takes the efficiency of absorption into account, 



EQE = 



no. of upconverted particles 
no. of incident photons 



^absXeqe: (3) 



where the constant n a bs = <7 a i> s (Ax) 2 , with sensitizer ab- 
sorption cross section o~ a bs, gives the number of created 
particles per photon incident in an area (Ax) 2 around 
a sensitizer molecule, and the efficiency measure Xeqe is 
defined below. Since the sensitizers cover only a fraction 
1/(1 + 1) of the total sample, we obtain: 

no. of upconverted particles 

EQE OC Xeqe = ^ J T , A ~ (4) 



T c t L/Ax 



X, 



qe i + r 



(5) 



where the second equality strictly speaking only holds 
for low densities. This assumption is fulfilled within the 
parameter space of consideration, ie. J ^> T c ,Td, as can 
be seen in [Figure 6\ where the density is on the order of 
10~ 4 <C 1. The external quantum efficiency x = Xeqe is 
used from here on as efficiency measure. 



III. NUMERICAL RESULTS 

In our numerical studies, we simulate a chain of sen- 
sitizer and emitter molecules on a lattice structure. The 
dynamic processes in the simulation are restricted to 
those shown in |Figurc 2} Thus, dynamics in the sublat- 
tice of sensitizers is simulated exactly, whereas dynamics 
of emitter molecules is subject to the constraints intro- 
duced above: an infinite triplet lifetime and no energy 
transfer from emitter to sensitizer. Due to these approx- 
imations, the kinetics only depends on the number of 
neighbouring sensitizers enclosed by emitters, but not on 
the number of neighbouring emitters. Therefore, we do 
not allow emitters to be placed next to each other on 
the lattice. To ease the analysis of numerical results, 
we study a chain of sensitizers where we place emitters 
equally spaced. The nonlinear two particle interaction 
coagulation is implemented by limiting the number of 
triplets to one per lattice site. Each triplet that reaches 



an emitter molecule, where it will stay until another 
triplet arrives and upconversion occurs, contributes to 
the upconverted fraction. 
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FIG. 3. Numerically obtained values for the external quan- 
tum efficiency, as function of emitter spacing 1, for the cre- 
ation rate T c = 10~ B (in units of the hopping rate, which 
is set to J = Is -1 throughout the paper, see Q). Each 
curve corresponds to a decay rate in the interval = 10 fc , 
k = — 1, . . . , —8. For decay rates with k — —5, . . . , —8, the 
curves almost coincide. Each curve exhibits a critical value 
Z c , where the external quantum efficiency is maximized. 

In all simulations, the rate of triplet hopping to the 
left and to the right was set to J — Is -1 . The actual en- 
ergy transfer rates, as given by the corresponding Miller- 
Abrahams rate, can be mapped to J = Is -1 by scaling 
all rates in our model (r c , Td) with the same factor, since 
the final results only depend on the ratio J/T c and J /Td- 
The scaling is described by the following equation, where 
Jtruei Ydjrue, Y ctrue refer to actual experimental values 
and r is the scaling constant: 



d,true 

T d 



- c,true 



(6) 



According to Figure 3| the efficiency can vary strongly 
with the emitter spacing I, which controls the sensitizer- 
to-emitter ratio, see ([!]), and is easy to control in experi- 
ments. It is therefore an ideal candidate for optimization. 
For large decay rates (Td — 10 k s~ 1 , with k = —1, —2, —3 



in Figure 3 ) , the efficiency has a pronounced peak at a 
very well defined value of I, which we name l c . Decreasing 
the decay rates increases l c . Finally, for very low decay 



rates, ie. k > 4 in Figure 3| the value of I at which the 
efficiency is maximal does no longer change. l c saturates 
due to particle loss by coagulation, and the efficiency is 
near-optimal in a large interval of I. This can be observed 
in Figure 4 where we plot l c as a function of the decay 



rate r^. It becomes apparent that the behavior of l c can 
be classified into two different regimes: on the one hand, 
for large T c and small Td, l c is constant as a function of 
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FIG. 4. Numerical results for the critical value l c of the 
emitter spacing (ie. the value for which x is maximal), as 
function of decay rate T^, for creation rates F c — 10 fc , k = 
— 1, . . . , — 8. The connecting lines are merely a visual aid. 
The two limiting regimes are clearly distinguishable: In the 
coagulation limited case, l c is constant, whereas the slope of 
l c is similar for all values of T c , in the decay limited case, and 
coincides with the analytical prediction (dashed line) given 
by p|. 



Yd, and mainly depends on the creation rate r c . In this 
regime, the efficiency is limited by coagulation, as shown 
below. On the other hand, for small r c and large T^, 
l c is independent of T c and scales according to a power 
law with Td, consistent with (12) derived below. In this 



regime, the efficiency is limited by decay. The transition 
between the two regimes occurs at the point in parame- 
ter space where the timescales of coagulation and decay 
are of the same order. We discuss this in section 5. 

In |Figurc 5 we plot the maximal upconversion effi- 
ciency xQc)- In order to highlight the regime of high effi- 
ciencies , we plot 1 — x(' c ) on a logarithmic scale. Again, 
we clearly distinguish between two different regimes: ei- 
ther the efficiency is constant as a function of and thus 
purely determined by T c (coagulation limited regime), 
or the efficiency is mainly determined by and ap- 
proximately given by (13) derived below (decay limited 
regime) . 



IV. ANALYTICAL RESULTS 

In the continuous limit I 3> 1 (i.e. spacing Ax between 
sensitizers much smaller than spacing I Ax between emit- 
ters), our model can approxi mate ly be described by the 
following mean field equatio n 1 31 * 41 ! 



FIG. 5. Maximal value of external quantum efficiency x(^c) 
for optimal emitter spacing, as function of decay rate Td, for 
r c = 10 fc , k — — 1, . . . , —8. The connecting lines are merely 
a visual aid. For small Td and constant r c , the efficiency is 
purely coagulation limited. For large F d it approaches the 
analytical prediction (dashed line) given by (13 1. 



where 4> is the density of triplets in the sensitizer, such 
that 4>dx is the probability to find a triplet in dx. The 
first term on the right hand side describes creation of 
triplets. T c describes the rate of triplet creation at a 
single molecule, where the molecule spacing is Ax. In 
the continuous system the creation rate is thus T c /Ax. 
For simplicity, we set Aa; = 1 (which is possible since 
the final result for the efficiency does not depend on the 
choice of Ax). The factor 1 — cf> accounts for the fact 
that only a sensitizer in its ground state can absorb a 
photon, which subsequently results in a triplet, whereas 
a sensitizer that is already in its T± state can not ab- 
sorb another photon. In principle, if the sensitizer has 
appropriate energy levels, one has to allow for photon 
absorption from the excited T\ state. However, the ab- 
sorbed energy is immediately lost, due to coagulation, 
and is therefore irrelevant for the kinetics. The second 
term describes triplet decay, e.g. via phosphorescence or 
non-radiative decay channels. The third term represents 
triplet diffusion with diffusion constant D = J(Ax) 2 = J 
for Ax = 1. Finally, the last term accounts for coagu- 
lation. ([7]) is almost identical to the upconversion rate 
equation that is commonly used to describe triplet-triplet 
annihilatiorP^. The only difference comes from the term 
cubic in the density, which correctly describes annihi- 
lation in one dimensio n 1 41 ! 42 !. The absorbing boundary 
conditions translate to 



<j>(x e ) = 0, 



(8) 



at 



TJl 



~Q^2 ~ Ctcoa 9 ( t ) ' (') 



where x e is the spatial coordinate of any emitter. 

A closed analytical solution to the nonlinear differen- 
tial equation Q is not available. However, in the low 
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is given by J times the first derivative of 4>, hence 



Xiqe 
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(11) 



FIG. 



Spatial dependence of mean triplet density for cre- 



where the factor 2 takes the two boundaries into ac- 
count, with a as above. The integrand corresponds to 
particle creation, see ([7]), and the limits 1/2 and I + 1/2 
of integration, with corresponding evaluation of the flux 
at x = 1/2, ensure that we count exactly I sensitizer 
molecules, and not I + 1. Due to the neglect of coagu- 
lation, the above treatment should correctly reproduce 
only the behaviour in the decay limited regime. To check 
this, we evaluate the maxima of the external quantum ef- 
ficiency x = l/(l+ l)Xig e , with Xiqe given by under 
the assumption 3> T c . Under the additional assump- 
tion <C J (i.e., that triplet hopping is much faster than 



ation rate r c 



10" 



and emitter spacing I = 30. Solid decay), we obtain the following analytical expression 



lines represent analytical results, whereas numerical data are 
indicated by symbols. The decay rate Yd = 10 fc s _1 is changed 
for each curve, with k = — 1, . . . , —6. The undermost curve 
corresponds to k = — 1, and the topmost to k = —6. 



density limit, <f) <C 1 or IT C <C T d + J/l 2 , the nonlinear 
term can be neglected, resulting in the linear differential 
equation 



6J\ 



1/3 



(12) 



dt 



r c (i 



dx 2 



(9) 



for the optimal emitter spacing in the decay limited 
regime. Using (11) and |5|, it is easy to verify that: 

As expected, both expressions agree with the numerical 
data shown in Figure 4 and Figure 5 (dashed lines). 



For convenience, we consider only positions between the 
first two emitters, i.e. < x < (I + l)Ax, since the solu- 
tion of ([9]) will be identical between all other emitters. In 
the steady state, ^ can be solved exactly for boundary 
condition ([8]), resulting in 



cosh [a(l+l- 2x/Ax)} 
cosh [a(l + 1)] 



(10) 



Figure 6 



this result is 
from the exact 



with a = y/ (r c + T d )/J. In 
compared with numerical data obtained 
many-particle simulation as described in |scction II| Not 
surprisingly, the prediction of the analytical equation is 
correct in the low density regime <f» <^L 1, which is realized 



for all curves shown in Figure 6 Remarkably, ( 10 ) also 
shows a good agreement with nume rics for r c = 10 i , 



T d = 10 s , which, according to Figure 4 is in the co- 



agulation limited regime. This shows that, also at small 
densities, the efficiency may be limited by coagulation - 
provided that the decay rate is even smaller. In this case, 
the efficiency is close to 1, see |Figure 5| 

The flux through the boundaries and the number of 
created particles define the internal quantum efficiency. 
According to Fick's law, the flux through the boundary 



V. SCALING ANALYSIS 

Above, we already derived the scaling behaviour in the 
low density limit. Since there is no analytical steady state 
solution of Q for the high density limit available, we 
present a more intuitive argument in this section, which 
works in both, low and high density regimes. 

Besides triplet escape to the emitter and subsequent 
upconversion, there are two other triplet loss mechanisms 
in the system: Decay, which governs the low density 
regime, and coagulation, which is predominant in the 
nonlinear high density regime, see |Figure 4"] According 
to the diffusion equation, a triplet placed in the center 
of the system needs the time r e = l 2 /J till it reaches 
the boundary. This defines the timescale of escape. The 
decay timescale is r d = 1/T d . Coagulation occurs on 
the timescale in which a triplet is created in the system 
with length I, thus r co = l/(Zr c ). If a particle exists al- 
ready in the system, there will be two particles after a 
time t co on average, and, consequently, coagulation can 
occur. The transition between decay and annihilation 
limited regimes is expected to take place when the re- 
spective timescales are of the same order: 



T d 



Td~ 



(14) 
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Clearly, the efficiency must depend on how the three 
relevant timescales mentioned above compare to each 
other. On the other hand, we have seen that the only 
relevant parameters that define the optimal structure are 
the creation rate in the coagulation limited regime, and 
the decay rate in the decay limited regime. Thus, in each 
regime, the fraction of triplets which escape and are up- 
converted is determined by only two timescales, which 
are r e and either or t co . Therefore, in the steady state 
the effective particle creation rate equals the sum of the 
two relevant loss rates: 



r c (i -</>) = r e + r. 



d,co • 



(15) 



The internal quantum efficiency is identical to the quo- 
tient of the escape rate and the effective rate with which 
particles are created. This implies: 



Xi- 



(16) 



In order to obtain an efficient system we demand r e < 
Td t co- Expanding in the small parameter up to first order, 
we find 



I 



X = 



l + l 



1 - 



Td.s 



(17) 



As we will show below, this intuitive equation yields the 
correct scaling behaviour for optimal emitter spacing and 
optimal external quantum efficiency, when analysed in 
the respective regimes. By evaluating the maxima of 
(17), we obtain the optimal emitter spacing l c . Solution 



of 



dx 
dl 



= 



(18) 



for I, evaluated in the limit of r c + Td -C J, yields l c and 
its scaling behaviour. 

In the decay limited low density regime we find 



1/3 



(19) 



for the scaling of l c , and by substitution in ( 17) we obtain 
the scaling law for the efficiency: 



i-x(4 d) )c< 



1/3 



(20) 



The scaling laws (19) and (20) agree with the analytical 
results (12) and (13). However, it is important to note 
that, while (17) gives the correct scaling behaviour, it 
can not be used to quantitatively predict absolute val- 
ues. The reason is the following: For evaluation of (17) 



we insert an escape rate, which is only an upper bound. 
In fact, the escape rate depends on the spatial coordinate, 
since a particle sitting close to the boundary will escape 
faster than a particle sitting in the center of sensitizer 
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FIG. 7. Optimal value of scaled emitter spacing yi = I c tV A 
as function of the scaled variable xi , see equations ( 23 1 and 
(24), for different values of r c (labels as in Figure 4 1. Due 
to the rescaling, the data for different r c collapse onto a 
single curve. Therefore, the scaled optimal emitter spacing 
is completely characterized by the single variable Xi. The 
dashed/dotted lines mark the behaviour y\ co ^ — 6 1//3 /xi and 
y ( ) _ gi/4 - n coa g U i a tion and decay limited regime, re- 
spectively. The dash-dotted line separates decay and coagu- 
lation limited regimes in parameter space, see (141 and (25 1. 



molecules. Averaging over all spatial coordinates would, 
in principle, yield the correct value. However, for aver- 
aging one has to know the density in order to give each 
contribution the correct weight, and the density in the 
coagulation limited regime is not available analytically. 
Thus, the absolute value of the inserted escape rate is 
not correct, but the scaling of l c with the system param- 
eters is, as it yields the correct exponent and agrees with 
the numerical data. 

Similarly, we obtain scaling laws in the coagulation 
limited high density regime, where an exact analytical 
result is not available. We find 



Kco) 



1/4 



(21) 



for the scaling of the optimal emitter spacing, whereas 
the maximal efficiency scales as 



x(ii co) ) « 



1/4 



(22) 



As above, we do not inherit the correct prefactors from 
(17). In order to determine these factors, we define 



rescaled variables 



Xi 



-.1/3 

_d 

,1/4 



(23) 
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FIG. 8. Scaled maximal efficiency y x , see (29 1, as a function 
of the scaled variable X ~y - SCO (|28n . Labels encode different 



see ([281. 

values of r c , as in |Figure 4| T he dashed/dotted lines mark 
the scaling prediction (|31[)/i 30 1 in the coagulation/decay lim- 
ited regime, respectively. The numerical data for different Y c 
and Yd collapse onto a single curve, except for very large Y c or 
Yd, where the assumption Y c + Yd <C J fails, ie. x x < 10 -0 ' 1 
and x x > 10 07 . The transition between decay and coagula- 
tion limited regime of the optimal efficiency occurs at x x = 1, 
where x x < 1 corresponds to decay limited optimal efficien- 



missing prefactor in (21 1 as 



z (co) 



8.7 

r7 



1/4 



(26) 



Concerning the coagulation limited efficiency, see (22 1, a 



fit of the numerical data shown in Figure 5 (for Yd — » 0) 
yields the prefactor: 



1/4 



(27) 



Correspondingly, if wc define scaled variables for \ as: 



V16J 




ior e 

33J 



1/4N 



and 



MS) 



1^ 



> (28) 



(29) 



the rescaled maximal efficiency as given by (13) and (27) 
is found as 



y« = 1 

in the decay limited regime, whereas 

4 CO) = h 



(30) 



(31) 



for the rates Y r and T r < 



in the coagulation limited regime. This is confirmed by 



m 



i r 1/4 



(24) 



for the optimal emitter spacing. According to equations 
(12) and (21), the rescaled emitter spacing yi as a func- 



tion of xi is constant in the coagulation limited regime, 
and given by y^' = (6J) 1 / 3 / 'x\ in the decay limited 



regime. Moreover, according to (14), the border separat- 
ing the decay limited regime from the coagulation limited 
regime turns into 



yi 



(25) 



i 1/3 

such that xi < y/ corresponds to the coagulation lim- 

1/3 

ited regime, and xi > y { to the decay limited regime. 
These predictions are confirmed by |Figure 7] where we 
show the same data as in Figure 4 (optimal emitter spac- 



ing as a function of F^) - but now expressed in terms of 
the rescaled variables xi and yi. Indeed, all data points in 
|Figurc 7] collapse onto a single universal curve, which, on 
the double logarithmic scale, is given by two straight lines 
corresponding to the two regimes discussed above. Fur- 
thermore, one of the lines agrees with the above predic- 
tion 2/t = 6 1 / 3 jx\ for the decay limited regime, whereas 

the value of the constant y^ = 8 1 / 4 ~ 1.7 observed for 
the coagulation limited regime allows to determine the 



Figure 8 where we plot the same data as in Figure 5 



in terms of the scaled variables x x and y x . A collapse 
onto a single universal curve can also be observed here. 
Note that some numerical values slightly deviate from the 
predicted curve, since these data points are not strictly 
in the limit + Y c <C J, for which the scaling equations 
were derived. 



VI. CONCLUSIONS 

Optimization of the sensitizer-to-emitter ratio can 
enhance the upconversion efficiency drastically. We iden- 
tified two characteristic kinetic regimes of the system, 
one where single particle decay limits the efficiency, and 
a second one where the nonlinear particle interaction 
coagulation (triplet-triplet annihilation) is the limiting 
factor. For each regime, we derived analytic expressions 
for the optimal ratio of sensitizers and emitters. The 
provided scaling laws for the optimal structure and 
efficiency reveal the universal behaviour of our model 
system. In future work, we would like to answer the 
question how the scaling changes in higher dimensions. 
Knowledge of the scaling laws of a 3D system, which 
most likely differ from the ID scaling laws shown here, 
will make it simpler to verify the suitability of a given 
material and optimize it for non-coherent upconversion. 
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